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By using the multiple scale method with the simultaneous introduction of 

a' 

multiple times, we study the propagation of long surface- waves in a shallow 
inviscid fluid. As a consequence of the requirements of scale invariance and 
absence of secular terms in each order of the perturbative expansion, we show 
that the Korteweg-de Vries hierarchy equations do appear in the description 
of such waves. Finally, we show that this procedure of eliminating secularities 
is closely related to the renormalization technique introduced by Kodama and 
Taniuti. 
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I. INTRODUCTION 



In 1967, Gardner, Greene, Kruskal and Miura Q, by making use of the ideas of direct 
and inverse scattering, showed that the Korteweg-de Vries (KdV) equation could be solved 
exactly as an initial value problem. Shortly after, Lax [Q] generalized these ideas to a general 
evolution equation of the type 

with K a nonlinear operator that could be written in the form 

K[u] = [B,L], (2) 

and with B and L self-adjoint linear operators depending on u. As a consequence of this 
formalism, Lax showed the existence of an infinite sequence of integrable partial differential 
equations of the form 

^ = K n [u] , n = l,2,..., (3) 

where 

K n =[B n ,L] (4) 



with 
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Bn = D 2n+1 + ]T (bjD 2j+1 + D 2j+1 bj^J (5) 



and 



L = D 2 + l -u. (6) 
o 

In these equations, bj are coefficients that can be determined 0, and D = d/dx. For n — 1, 
the KdV equation is obtained. For n > 2, we have the higher order KdV equations. This 
infinite sequence of integrable nonlinear partial differential equations form the so called KdV 
hierarchy Q. Along with KdV, they are all integrable by the inverse scattering method, 



and they have the same integrals of motion as KdV. However, in contrast to KdV, which 
has been shown to govern the nonlinear long-wave dynamics of general dispersive systems, 
the physical relevance of the higher order equations of the KdV hierarchy has remained, up 
to now, obscure. 

We will show in this paper that the equations of the KdV hierarchy do appear in the 
description of physical systems. More specifically, we will first argue that the same ampli- 
tude that satisfies the KdV equation, also satisfies the higher order equations of the KdV 
hierarchy, each one in a different time-scale. As a consequence, we will then show that 
this property prevents divergent solutions to the physical problem. The infinite sequence 
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of time-scales are defined from the basic time variable t by T\ = e?t, r 3 = eat, r 5 = e^t,..., 
where e is a parameter satisfying e<l. The primary reason for introducing them is that 
they may be used as a tool to eliminate the secularities that appear in the problem when a 
perturbative solution for the higher order terms of the amplitude is searched. However, when 
a function of a single time variable F(t) is extended to a function of several time variables 
F(tx, T3, T5, ..), the problem of knowing the allowed evolution in each one of the different 
time-scales is posed ||. Indeed, as we are going to see, this evolution is not arbritary. It 
is determined by a scale invar iance requirement, which must hold to ensure the ordering of 
the expansion, and by the compatibility condition 

d 2 F d 2 F 
dr 3 dr 2n+1 dT 2n+1 dT 3 ' 

If, for a physical system, the evolution of a certain amplitude is governed by the KdV 
equation in the time T3, we will show that the above compatibility condition constraints 
the evolution on higher order times of that amplitude to be governed by the higher order 
equations of the KdV hierarchy, leaving only one free parameter at each order which is 
related to the possibility of redifining the time r 2n +i by a multiplicative factor a 2n +i- By 
choosing specific values for the free parameters, the task of secularity elimination can then be 
accomplished. As a consequence, our construction of a bona fide perturbative expansion will 
provide a link between the equations of the KdV hierarchy and the evolution of a physical 
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quantity. The physical system studied here, that of long surface water- waves, is a kind of 
classical system where our ideas are well illustrated. However, the results obtained are, to a 
certain extend, model independent, and in this sense the study made here can be considered 
as representative of a general physical situation. 

The paper is organized as follows. In section II, we obtain the basic equations describing 
surface water waves in the multiple space and time formalism. In section III, the first three 
evolution equations are obtained. In section IV, by using the multiple time formalism, we 
examine how the symmetry of the time derivatives determines the evolution of the wave 
amplitude in any time scale. The use of the KdV hierarchy equations to eliminate the 
secularities in the evolution of the higher order terms of the wave amplitude is discussed in 
section V. In section VI, we consider the case of a traveling-wave solution to the equations of 
the KdV hierarchy, and we show that the method of eliminating secularities by using these 
equations can be related to the results of Kodama and Taniuti [|J, where a renormalization 
technique was introduced to obtain a secular free perturbative expansion. And finally, in 
section VII, we summarize and discuss the results obtained. 

II. THE MULTIPLE SPACE AND TIME FORMALISM FOR WATER WAVES 

We consider a two-dimensional inviscid incompressible fluid in a constant gravitational 
field. The space coordinates are denoted by (x, z) and the corresponding components of the 
velocity v by (u,w). The gravitational acceleration g is in the negative z direction. The 
equations describing such a fluid are 0: the incompressibility equation 



V -v = 0, 



(7) 



and the Euler equation 




(8) 



where p is the fluid density, p is the presure and k is an unit vector in the z direction. If we 



assume the fluid to be irrotational, 
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Vxw = 0, (9) 
and consequently a velocity potential <fi can be introduced through 

v = V(t>. (10) 

Using this definition, Eq.(^) becomes 

V 2 = O, (11) 
while Eq.(HD, after an integration, reads 

— = -Jt (t) - ~ 9Z > (12) 

with p an integration constant. 

We consider the case of a fluid of height h, limited above by a passive gas exerting a 
constant pressure po on it, and let the upper surface to be described by 

z = ((x,t). (13) 

The kinematic boundary condition at this surface is then written in the form 

DC 

■^ = Ct + = (14) 
with the subscripts denoting partial derivatives. However, there is also a dynamic boundary 



condition obtained from Eq.(12), which reads 

1 



>' + 2 



,) 2 + (<fc) 2 +gp = 0, (15) 



onz = ((x,t). Finally, the lower boundary is supposed to be a rigid horizontal flat botton, 
localized at z = —h. In this case, the corresponding boundary condition implies that the 
normal velocity of the fluid must vanish: 

2 = O. (16) 

We now wish to consider the long-wave in shallow-water approximation to the above 
equations. This may be done, for instance, by using the the reductive perturbation method 



of Taniuti 0, which introduces slow space and time variables. The slow space variable is 
given by: 

£' = e*z. (17) 
Next, we introduce a very general set of slow time variables J7|: 

n = eh , T 3 = eh , r 5 = eh . (18) 
In addition, we expand ( and in a suitable power series in the parameter e: 

C = 6C = e(Co + eCi + e 2 C 2 + ...) (19) 

cf) = e ^</> = (0 O + e 0i + e 2 2 + • • •) • (20) 

Introducing these expansions and the slow variables into the water wave equations, we 
obtain: 

e^> ee + (p zz = , -h< z < eC (21) 

Z = , z= -h (22) 

2 = eCn + e 2 Cr 3 + e 3 Cr 5 + • • • + e 2 0?'4' > * = < (23) 

2g( + 20 T1 + 2e0 T3 + 2e 2 T5 + . . . + e^'k + kk = , ^ = eC • (24) 

These are the basic equations we are going to use to obtain the evolution equations describing 
the system. 

III. THE FIRST EVOLUTION EQUATIONS 

Equation ( |2TJ ) for the velocity potential can be solved for the boundary condition (P2"|). 
independently of equations (|23D e (f24l). At order e°, the solution is 
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(25) 



with T = J^iC'i r i) r 3> r 5; • • •) an arbitrary function. At order e 1 , it is given by 

0i = - + M ^f'f + G , 

with = (/(£'; Tx, r 3 , r 5 , . . .) another arbitrary function. At order e 2 , we get 

2 = -— (z 4 + 4hz 3 - 8h 3 z) F m - - (z 2 + 2hz) g ee + ft , 

with n = H(%\ T ii r 3; r 5> • • •) another arbitrary function. At order e 3 , the solution is 

1 



(26) 



(27) 



z 6 + Qhz 5 - 40/i 3 z 3 + lQh e ) F m 



+ y (V + 4hz 3 - 8h 3 z) Q m -\{z 2 + 2hz) H ee + X, 



(28) 
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with X = Z(£'; r 1; r 3 , r 5 , . . .) again an arbitrary function. We could proceed further and 
calculate 04, 05, . . . . But, we stop here since this is all we are going to need. 

We pass now to eqs.(^) and fl2~4|), and solve the first few orders in our perturbative 
scheme. At order e°, Eq.(^) gives nothing, while Eq.(|24|) reads 



Substituting O , we get 



>0n 



1 







(29) 



At order e\ Eq.flU is 



Co — F- n • 

9 



F T1T1 - ghT^i = . 



(30) 



(31) 



Defining 



gh = c 2 



(32) 



with c a velocity, its solution can be written in the form 



f = F(Z'-cr 1 )+G(e + cT 1 ) 



where the sign — (+) refers to a wave moving to the right (left) with a velocity c. For 
definiteness, let us choose the wave moving to the right, and let us define a new coordinate 
system by 

£ = ? - cn . (33) 

Therefore, 

T TX = -cFz , (34) 



and from Eq . (|30|) we get 



Consequently, 



Co = ■ (35) 
c 



Con + cCo f = , (36) 



which is the usual linear wave equation. 

Now, using eqs.flTTD and (|T8|), we can rewrite Eq.(p3|) in the form: 



£ = (x ~ ct) , (37) 

from which we see that £ is a coordinate system moving with a velocity c of the linear waves, 
in relation to the laboratory coordinate x. Therefore, from now on, all dependent variables 
will be assumed to depend on (£', n) through £ only, which automatically implies that they 
satisfiy Eq.(|3J). 

At order e 1 , Eq.(|2^) gives 

2g( 1 - 2c£f + 2 JV 3 + = , (38) 

where we have already used Eq.(|34]) for Q. Derivating in relation to £, and using Eq. (p5|) , 
we obtain 

cCi« - hGtt = -Cor 3 - T^CoCo? • (39) 



Now we go to the order e 2 . From Eq.(^), we have 



c ch 2 

c&* - hQx = Co T3 + 2-CoCo § + — Co« € • (40) 
Substituting into Eq.(^), we get 

3c ch 2 

Co7 S + 2^CoCoe + -g-Coece = o, (41) 

which is the KdV equation. It is the first nontrivial equation in the KdV hierarchy. At 
this point we can see the importance of introducing multiple times: while the linear wave 
equation fl56] ) describes the evolution of Co i n the time r 1; the KdV equation describes the 
evolution of (q in the time r 3 . In other words, the phenomena described by each one of these 
equations occur in different time scales. And this will be the case for all the higher order 
equations of the KdV hierarchy. Now, from Eq . (|24[) at order e 2 , we obtain 

2g( 2 + 2c 2 CoCo« - 2cft 5 + 2Q TZ + 2T Th + + c 2 (o^ = . (42) 

Derivating in relation to £, and using Eq . fl35|) , it becomes 

h 

<2§ - hH^ + -Qfa + (Coe^ + CoGix) = -Cor 5 - 2c/iCo^Co?c - c/iCofe • (43) 



We pass now to the order e 3 . From Eq . (|23"D , we have 
< 2 £ - kHz - (Co% + Cof^e) 



c h 3 2 

-r (C05C1 + C0C15) - = Cor 5 + Y^ c/l4 Co(50 • (44) 

Equations (|^) and fl44|) can be combined to yield an evolution equation involving (j T3 , ( 0T5 
and ^ T3 : 

2 4 

2C0r 5 + J^c/l Co(5?) + 2cftCoeCo« + c/l CoCoC« + Clr 3 

+ £ (Co € Ci + CoCi € ) + -Qir, + 2 (Co% + Co^) + y £(40 = . (45) 
Now, making use of Eq.([y]) to describe T T3 , Eq.(|38|) can be put in the form 

f'-lft^' + T^- (46) 



We can then use this equation to eliminate Q% from Eg. (|45|) . The result is 



where 



S(Co) = -Cor 5 - ^ch*( om - ^ch( ( om - ^chC^Cm + ^Co 2 Co€ • (48) 



Equation (|47|), as it stands, can not be viewed as an evolution equation for d in the time T3 
because the non-homogeneous term involves (o T5 , and the evolution of (0 in r 5 is not known 
up to this point. Moreover, the term proportional to Co(5£) m 5"(Co) is a resonant term, that 
is, it is a secular producing term to the solution (1 M. We will show in the next sections 
how the evolution of Co m r 5 may be used to cancel the secular term. However, let us first 
state a mathematical result related to the KdV hierarchy, which will be used later. 

IV. THE SYMMETRY OF TIME DERIVATIVES AND THE 
KORTEWEG DE VRIES HIERARCHY 

As is widely known M, the KdV equation 

Ct = 6CC« - C«« , (49) 

is a Hamiltonian system, and can be written in the form 

d (8R 



with the Hamiltonian H given by 



H 



-co 



(Cx) 2 _j_ ^3 



dx. (51) 



2 

On the other hand, as is also well known, the KdV equation has an infinite sequence of 



conservation laws |]|TI|, with the integrals of motion given by 



+00 



J n [C] = / T n dx , 71 = 0,1,2,3,... , (52) 
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where T n , the conserved density, is a polynomial in CXx, Cxx, etc. Their explicit form can 



be obtained through a perturbative scheme developed by Miura, Gardner and Kruskal [jTO 
which is based on expansions on a small parameter e. In this method, each integral of motion 
I n appears in a different order of the perturbation parameter e. Now, the infinite sequence 
of equations of the KdV hierarchy can be expressed in terms of these integrals of motion 
according to 



For n = 1, the integral I\ coincides with the Hamiltonian given by Eq.([5l|), and the KdV 
equation is obtained, which is the first nontrivial member of the KdV hierarchy. From the 
explicit form of /„ [1C], it is easy to show that each one of Eqs.(|53D is invariant, up to a 



Galilean term, under the transformations 

£ = e 5(x-ct) (54) 

T 2n+1 =e n+l *t ; n = 1,2,3,... (55) 



C = eCo- (56) 
In these new coordinates, called slow variables, Eq.(^) is written in the form 

Co — = m (ji) ■ (57) 

Notice that an infinite sequence of slow time variables r 3 , r 5 , r 7 , ... was introduced in the 
above transformations. As a consequence, any dependent variable F(x; t) will be a function 
of these multiple times: 

F(x',t) = F(x',T 1 ,T 3 ,T 5 ,...). (58) 
Therefore, the time derivative of F(x; t) turns out to be 

8F ( 3 8 5 d 7 d \ „ 

a = r^ +£ ^ +e, ^ + -r- (59) 

n 



Let us now return to our physical system. As we have already seen, the evolution of Co 
in the time r 3 was found to be the KdV equation 

C0r 3 = ttlCo«€ + Mofa , (60) 

with 

ch 2 3c , . 

a, = — , ^ = - h . (61) 

Let us now examine how the evolution of Co in times t^,,tj,... can be obtained. The first 
crucial point of this analysis is to note that, for the evolution of Co in the time t$, the two 
terms in the r.h.s. of Eg . (|60"D exhaust all possible terms such that, when transforming back 
from the slow to the laboratory variables, the resulting eguation does not depend on e. 

In the very same way, the evolution eguation of Co i n the time T5 must only involve terms 
presenting the same property as the terms in the KdV eguation. If this where not so, the 
ordering of the perturbative series would not be ensured. In other words, this evolution 
eguation must be formally the same in the slow (£, T5, Co) as well as in the laboratory coor- 
dinates (x, t, C), establishing thus a scale invariance reguirement. A simple analysis shows 
that the only possible terms are: 

C0r 5 = «2C0(5O + /?2CoCo5« + (#2 + I2) CofCotf + ^Co^O? , (62) 

where <y.i,fii, 72,^2 are constants. Now comes the second crucial point: the coefficients 
&2, P2, 72, #2 are not completely arbitrary if Co is to satisfy also the KdV eguation in the time 
r 3 . Instead, they are constrainted by the relations arising when the natural (in the multiple 
time formalism) compatibility condition 

(CorX = (ConsXs (63) 

is imposed. This condition only makes sense in the multiple time formalism, since otherwise 
it would be redundant and would lead to a trivial identity. A tedious but straightforward 
calculation shows that the relations arising from this condition are: 
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«2 3 «i ' «2 3 «i ' «2 6 
Substituting into Eq.(B^) leads to 



(64) 



Cor 5 = «2 



Co(5€) 



I (f ) CoCo ^ + y (^) + 1 (^) 



(65) 



which is the 5 i?t order equation of the KdV hierarchy. Different choices of the free parameter 
a 2 correspond to different normalizations of the time T5. In particular, for a 2 = 6 it acquires 



the canonical form 12 



Cor 5 = 6Co(5§) + 10aCoCo$« + 20aCo?Co« + 5a C0C05 , (66) 

with a = fix/ ai. 

This procedure can be extended to any higher order. In other words, if (q satisfies the 
KdV equation in the time T3, by making use of the scale invariance described above, we can 
first select all possible terms appearing in the evolution of Co m the time r 2n+ i. Then, by 
imposing the general time symmetry condition 

d 2 Co = <9 2 Co 
dT 3 dr 2n+1 dT 2n +idT 3 ' 

the parameters a n , f3 n , j n , . . . appearing in that equation can be uniquely determined in 
terms of a\ and fii, and the resulting equation is found to be the (2n + l) th order equation of 
the KdV hierarchy. As before, a n is left as a free parameter responsible for the arbitrariness 
in the scale of the time r 2n+ i. One gets easily convinced of the validity of this result by 
performing the calculation explicitly for the first few equations, but we do not intend here 
to give a general proof valid for any order. 

Before proceeding further, let us remark that one could a priori envisage, in place of 
Eq.(62) a new one including also a dependence on (1, and still keep the agreement with the 
scale invariance requirement. However, when introducing this back into the physical equation 
(60), the compatibility condition (63) makes all the supplementary terms to vanish. 
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V. HIGHER ORDER EVOLUTION EQUATIONS: THE ELIMINATION OF THE 

SECULAR PRODUCING TERMS 



Let us now return to the evolution equation for (j in the time T3, which is given by 
Eqs.(|4"7|) and (fjgp. As we have seen, there were two problems related to it: the evolution of 
Co in r 5 , and the secular producing term. The last section gave the clue to the first problem. 
As we have seen, the evolution of Co in r 5 can not be chosen arbitrarily, but it is restricted 
to the KdV hierarchy equations. We show now how the secular producing term Co(5£) can be 
eliminated. A comparison between Eq. (|48|) for S(( ) and the 5 th order KdV equation (|65|). 
immediately shows that the choice 

1 Q 

a 2 = -—ch A (68) 
360 y ' 

has the property of eliminating the terms 

19 

± Cor 5 - ^ c/l4 Co(50 (69) 

from the non-homogeneous term S'(Co)- Consequently, S^Co) acquires the form 

9 5 63 c 

S = ^ ch ^Com + g c/i Co f Co« + y^Co 2 Co € , (70) 

and we see that it does not present a secular producing term anymore. 

It should be noted that the resonant term is the linear term of S((o). This property 
holds in any higher order of the perturbative scheme Q]. Therefore, the elimination of the 
secular producing term by choosing an appropriate value for t he free parameter a n in the 
higher-order equation of the KdV hierarchy, also remains possible in any higher order. 

VI. TRAVELING WAVE SOLUTIONS: THE RENORMALIZATION OF THE 

SOLITON VELOCITY 

We have already assumed, in the context of the multiple time formalism, that 

Co = Co(£;t 3 , t 5 ,t 7 , . . .) , (71) 

14 



where 

e = e'-cri, (72) 

so that Co automatically satisfies the linear wave equation 

Con + cCo 5 = . (73) 

Our concern now will be the solutions to the higher order equations of the KdV hierarchy. 
Hereafter, for simplicity, we will assume that c and h can be set equal to unity, so that the 
KdV equation reads 

Cor 3 = aiCoce? + ftCoCae > ( 74 ) 

with 

13 

«i = g , Pi = g ■ ( 75 ) 

We now look for a traveling wave solution Co that satisfies this equation in the time t$, but 
that satifies also the higher order equations of the KdV hierarchy in the times T5, T7, . . . . 
This multi-solution can be written in the form 

Co = y sech 2 (fcA) (76) 

with the argument A given by 

U2 1,4 u6 

1 J- ™ ™ /—.-\ 

A= f-sr > -ir t -^r Y — ■■ (77) 

and with Ai, A 2 , A3, . . . parameters depending respectively on the constants a%, «2, «3, 
This argument is to be interpreted in the following way: when the evolution equation under 
consideration is concerned to the time T2 n +i, the relevant argument is 

i=l A i 

with 9 n a phase involving all the remaining terms of the sum. Then, for the KdV equa- 
tion ([74]), A\ = (l/ a = 6, and the solution is 
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k 2 

Co = — sech 2 
so 3 



kU-jT 3 -6 3 (k;r 5 ,...) 



(78) 



where 83 is a phase. On the other hand, the travelling wave solution satisfying simultaneously 
the KdV equation (|74j) in the time T3, and the 5 th order KdV equation fl65f) in the time T5, 
is given by 



k 2 

Co = ysech 2 



kU- yr 3 - ^-r 5 - 9 5 (k; r 7 , . . .) 



(79) 



and so on, to any higher order. 

Once we have found the solutions to any member of the KdV hierarchy, let us return to 
the equation for Ci in the time t%, which is now (c = h = 1) written in the form 



Ci7 S + |(CoCi) € + ^Ci«f = 5'(Co), 



(80) 



with 



S((o 



9 
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24 -CoCoe« + gCoeCoK + j^Co Co« • 



(81) 



It is important to remember that the linear term of 5(Co), which is the secular producing term 
in the equation for Ci, was eliminated through the use of the 5 th order KdV equation (|65D, 
with «2 given by Eq. (|68|) . Therefore, the Co appearing in Eqs.(|80|-pTD is that given by Eq. (^9|) , 
since it must be a solution of KdV as well as of the KdV hierarchy 5 th order equation fl65|). 
The secular-free solution Ci can then be found by solving the linear equation (|80|-|3TD. This 
solution can be found in Ref. ||. 

Finally, let us take the solutions Co? an d let us make the transformation back from the 
slow (£, r, Co) to the laboratory coordinates (x, t, (). For the case of solution (]78D to the KdV 
equation, we get 

k 2 



c 



-e sech 



kt 1 ' 2 (x-V 3 t) 



(82) 



where 



V* 



k 2 



(83) 
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is the solitary-wave velocity in the laboratory coordinates. For the case of solution ( f79|) to 
the 5 th order KdV equation, we get 

C = %-e sech 2 [jfc e 1/2 (x - V 5 tj\ , (84) 

where now 

V, = V 3 + e 2 L . (85) 

Since A2 depends on the parameter 012, to choose «2 means to define the velocity renor- 
malization itself. In this sense we can say that there is a unique velocity renormalization 
leading to a secular-free perturbation scheme for (\. From these properties, we can see now 
that this method is equivalent to the renormalization technique developed by Kodama and 
Taniuti ||, in which the secular-free higher order effects were also given by the renormal- 
ization of the KdV soliton velocity. Moreover, in the same way as in the method of Kodama 
and Taniuti, if higher order scales are introduced, it is possible to continue the secular-free 
perturbation to higher orders by using the higher order equations of the KdV hierarchy. And 
in general, depending on the order considered in the perturbative scheme, the renormalized 
solitary-wave velocity in the laboratory coordinates is given by 

n k2« 

VW = c + £e 4 — • (86) 
The analysis for multi-soliton solutions proceeds in a similar fashion as well. 



VII. FINAL REMARKS 

By using the reductive perturbation method of Taniuti ||, with the introduction of 
an infinite sequence of slow time variables r±, r 3 , r 5 , . . . , we studied the propagation of 
long surface waves in a shallow inviscid fluid. The three main ingredients of our analysis 
were: (i) the scale-invariance argument, which restricts the possible forms of the evolution 
equations, and which is necessary for the coherence of the perturbative expansion; (ii) the 
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compatibility condition (63), which appears when extending a function of a single time- 
variable to a multiple time-variable; (iii) the secularity elimination procedure, without which 
the perturbative expansion would be meaningless. By using (i) and (ii) we have shown that, 
if the amplitude Co satisfies the KdV equation, it also satisfies the higher-order equations of 
the KdV hierarchy. Then, by using the 5 th order equation of the hierarchy with a properly 
chosen at2, we have shown that the secular producing terms of the equation for d in the 
time r 3 could be eliminated, a result which can be extended up to any higher order. As 
the secularity elimination is mandatory for a physical theory, being essencially a finiteness 
requirement, the results obtained in this paper allowed us to give a physical meaning to 
the KdV hierarchy equations. Thereafter, by considering a solitary wave solution, we have 
shown that the elimination of secularities through the use of higher order KdV equations 
corresponds, in the laboratory coordinates, to a renormalization of the soliton velocities, as 
obtained previously by Kodama and Taniuti ||. 

The study of long-waves in shallow water is representative of a wide class, that of the 
weak nonlinear, dispersive systems, where the KdV equation has a kind of universal char- 
acter. In this sense, we can say that the results of the present paper do not depend on the 
specific physical system under consideration, or, which is the same, on the specific form of 
the basic equations, except for the values of the coefficients in the perturbative expansion. It 
is, therefore, legitimate to conjecture that they might be extended to the above mentioned 
larger class of systems. 
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